Mid-rapidity charge distribution in peripheral heavy ion collisions 
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The charge density distribution with respect to the velocity of matter produced in peripheral heavy 
ion reactions around Fermi energy is investigated. The experimental finding of enhancement of mid- 
rapidity matter shows the necessity to include correlations beyond BUU which was performed in the 
framework of nonlocal kinetic theory. Different theoretical improvements are discussed. While the 
in-medium cross section changes the number of collisions, it leaves the transferred energy almost 
unchanged. In contrast the nonlocal scenario changes the energy transferred during collisions and 
leads to an enhancement of mid-rapidity matter. The renormalisation of quasiparticle energies can 
be included in nonlocal scenarios and leads to a further enhancement of mid-rapidity matter distri- 
bution. This renormalisation is accompanied by a dynamical softening of the equation of state seen 
in longer oscillation periods of the excited compressional collective mode. We propose to include 
quasiparticle renormalisation by using the Pauli-rejected collisions which circumvent the problem 
of backflows in Landau theory. Using the maximum relative velocity of projectile and target like 
fragments we associate experimental events with impact parameters of the simulations. For periph- 
eral collisions we find a reasonable agreement between experiment and theory. For more central 
collisions, the velocity damping is higher in one-body simulations than observed experimentally, 
because of missing cluster formations in the kinetic theory used. 
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I. INTRODUCTION 

Numerical simulations based either on the Boltzmann 
equation (including the Pauli blocking it is often called 
the BUU equation Q) or on the closely related method 
of quantum molecular dynamics (QMD) 0,0] are exten- 
sively used to interpret experimental data from heavy ion 
reactions. Due to their quasi-classical character, they 
offer a transparent picture of the internal dynamics of 
reactions and allow one to link observed particle spectra 
with individual stages of reactions. Naturally, simulation 
results are only approximate. For example, BUU simula- 
tions fail to describe the energy and angular distribution 
of neutrons and protons in low and mid energy domain 



In particular, the formation of a neck-like structure 
in peripheral heavy ion reactions and the impact on 
the fragmentation mechanism and production of light 
charged particles has been discussed for a couple of years 
0-0. It has been suggested that this neck instability 
can be important for the fast decomposition of matter 
and is probably neutron rich |2Cf| . Theoretical inves- 
tigations suggest that the neck is not formed in usual 
heavy ion simulations starting from the Landau equation 
pl|-p3| or BUU equations p4|, p5| in cluding additional 
mean field fluctuations derived in ]2fj|]27| and tested [S . 
The inclusion of fluctuations in the Boltzmann (BUU) 
equation has been investigated in 
Boltzmann-Langevin pictures |3l] |3C 

We will take here the point of view that the fluctua- 
tions should arise by themselves in a proper kinetic de- 




resulting in 



scription where all relevant correlations are included in 
the collision integral. The collision will then cause both a 
dephasing and fluctuation by itself. This procedure with- 
out additional assumption about fluctuations has been 
given by the nonlocal kinetic theory ]37|-|39[| and applied 
to heavy ion collisions in |^ . We claim that the de- 
rived nonlocal off-set in the collision procedure induces 
fluctuation in the density and consequently in the mean- 
field which are similar to the one assumed ad-hoc in the 
approaches above. As a first step to verify this we will 
investigate the formation of a neck and the mid-rapidity 
emission pattern. Indeed as we will demonstrate, the 
neck is much more pronounced if a nonlocal kinetic the- 
ory is used. 

Recent INDRA observation shows an enhancement of 
emitted matter in the region of almost zero relative ve- 
locity which means that matter is stopped during the 
reaction and stays almost at rest [ ^9p^ ]. This enhance- 
ment of mid-rapidity distribution can possibly be asso- 
ciated with a pronounced neck formation of matter. A 
pioneering work on describing such mid-rapidity emis- 
sion has been done by E. Galichet and F. Gulminelli 
The main problem is to find a proper selection of exper- 
imental data fisjl . The comparison between data and 
simulation has been performed according to cuts in the 
transverse energy. This was possible since the applied 
coalescence model allows to get rid of the Fermi motion. 
Since we want to omit coalescence we select the proper 
comparison with respect to the maximum velocity of the 
projectiles with respect to the ratio of transverse to to- 
tal energy. The latter ratio gives a nearly Fermi motion 
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independent scaling. 

We want to investigate here the peripheral heavy ion 
collisions and want to discuss different theoretical im- 
provements of the BUU simulations. We start from the 
nonlocal BUU equation which includes from microscopic 
derivation the effect of binary correlations on the colli- 
sion process. By this way we obtain a nonlocal off-set of 
the collision partners which account for the readjustment 
of the trajectories according to the virial corrections, 
which would be for hard spheres the excluded volume. 
This off-set induces mean-field fluctuations similar to the 
above mentioned improvements of BUU. We will give a 
combined picture of nonlocal off-sets and quasiparticle 
renormalisation which leads to the consistent inclusion 
of binary correlation on the collision integral and mean- 
field fluctuations. Section II represents the formerly de- 
rived nonlocal shifts in terms of an intuitive picture and 
discusses symmetries and implementation in numerical 
codes. The quasiparticle renormalisation is suggested in 
a novel form using the excluded events by Pauli - block- 
ing. In section III we present the numerical results and 
compare different approximations with the experimental 
data. Finally in section IV we summarise and give some 
outlook. 



II. THEORETICAL PRELIMINARIES 

The need for nonlocal corrections can be stimulated 
by discussing the scattering of two particles as superpo- 
sitions of wave packets and similar used in . The 
asymptotic wave packet after scattering can be written 
for large distance x from scattering centre 
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with the scattering amplitude /(re, cos 0) where 9 is the 
angle between the relative momenta before, re, and after 
the collision, re/. We proceed now and expand the scat- 
tering amplitude around the final difference momenta re / 
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The derivative of the phase 8 leads now to the definition 
of the effective space shifts A? and the time shift A|| 
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where we denoted the shifts corresponding to the direc- 
tion of re as || and _L Rewriting (111) we obtain 



We observe three effects of scattering on the asymptotics: 
(i) a genuine time delay Aii, (ii) an effective displace- 
ment of the two colliding particle of with respect to 
the centre of mass and (iii) a modification of scattering 
probability f = T{1 + (re - K/)V K |/(p, cos 0) 

The effect of non-local collisions on the dynamics of 
heavy ion reactions has been studied already within a cas- 
cade model p8| . For a simple hard-sphere approximation 
of nucleon-nucleon collisions, Halbert has demonstrated 
that density patterns of 20 Ne+ 238 U reactions are sensi- 
tive to local or non-local treatment of collisions. Malfliet 
p9[ also found disturbing that all dynamical models rely 
more or less on the use of the local approximation of bi- 
nary collisions, because the local approximation neglects 
a contribution of the collision flux to both material re- 
lations which control the hydrodynamic motion during 
the reaction, the compressibility and the share viscos- 
ity. To include the collision flux, Malfliet incorporated 
hard-sphere non-local collisions into the BUU simulation 
code. Recently, this approximation has been used by Ko- 
rtemeyer, Daffin and Bauer |50| . 

The hard-sphere approximation of a non-local collision 
is sufficient for the above mentioned discussions of trends, 
however, it cannot be used in realistic studies. This ad 
hoc approximation has been used not only for its sim- 
plicity but also because of lack of a first principle theory 
offering quantum mechanical displacements which would 
generalise the classical hard-sphere displacements pro- 
posed by Enskog. As far as we know, till recently there 
was no non-local theory of binary collisions devoted to 
the nuclear matter. In literature, there are closely re- 
lated quantum theories of binary collisions developed for 
moderately dense gases |5l]]5^], these, however, treat non- 
local collisions via gradient contributions to the scatter- 
ing integral. The gradient form is suitable for hydrody- 
namic expansions studied in the chemical physics, but is 
very inconvenient for numerical simulation and thus have 
never been employed for heavy ion reactions. 

Recent theoretical studies have filled this gap in theory. 
Danielewicz and Pratt |53j pointed out that the collision 
delay can be used as a convenient tool to describe the 
virial corrections to the equation of state for the gas of 
quasiparticles. Although their discussion is limited to 
the equilibrium, it marks a way how to introduce virial 
corrections also to dynamical processes. The kinetic 
equation for quasiparticles with non-instantaneous and 
non-local scattering integral has been derived in J37 54 
as a systematic quasi-classical limit of non-equilibrium 
Green's functions in the Galitskii-Feynman approxima- 
tion. It has been shown that the gradient corrections to 
the scattering integral can be rearranged into a form of 
a collision delay and space displacements reminiscent of 
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classical hard spheres, i.e., into a form suitable for nu- 
merical simulations. 

In this contribution we will put these ideas of nonlocal- 
ities on the firm ground using the quantum kinetic equa- 
tion with nonlocal scattering integrals which was derived 
from quantum statistics J37],|54| to show how the effect 
of nonlocalities play a role in simulations of heavy ion 
reactions and compare them with experiment. 



A. Nonlocal kinetic theory 

The scattering integral of the non-local kinetic equa- 
tion derived in p7| corresponds to a following picture of 
a collision as seen in Figure |l} 
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FIG. 1. A nonlocal binary collision according to Eq. (^). 

Assume that two particles, a and 6, of initial momenta 
k and p start to collide at time instant t being at coordi- 
nates r a and r b . Due to a finite range of the interaction, 
at the beginning of collision particles are displaced by 
ri, — r a = A bc . The collision has a finite duration A t , 
i.e., it ends at t + A t . During the collision, both particles 
move so that their end coordinates differ from those at 
the beginning, r' a — r a = A a and r' b — rb = A&. The parti- 
cle a transfers a momentum q to the particle 6, therefore 
their relative momentum changes from k = — p) to 
k! — h(k—p)—q. Their sum momentum is modified by an 
external field acting on the colliding particles during the 
collision going from K = k + p to K' ~ k + p + Ak- The 
same field changes the sum energy of colliding particles 
from E = e a + e b to E' = e' a + e' b = e a + e b + A E - 

The values of A's are given by derivatives of the scat- 
tering phase shift <j) = Im lnTn(f2, k,p, q, t, r), 
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Note that energy Q enters as an independent quantity 
so that one needs to know the scattering phase shift out 
of the energy shell. The on-shell energy, fi = E, is sub- 
stituted after derivatives are taken. Since experiments 
provide us only with the on-shell values of the scattering 
phase shift </>, set (||) of A's cannot be derived directly 
from experimental values but only indirectly via off-shell 
T-matrix constructed from the interaction potential. Be- 
low we show how to circumvent this complication. 



B. Instantaneous approximation 

It is our intention to incorporate these features of col- 
lisions into the BUU (or QMD) simulation codes. The 
selfconsistent evaluation of all A's for all collisions would 
be too demanding. We employ two kinds of additional 
approximations. First we neglect the medium effect on 
binary collisions, i.e., use the well known free-space T- 
matrix. Second, we rearrange the scattering integral into 
an instantaneous but non-local form. This instantaneous 
form parallels hard-sphere-like collisions what allow us 
to employ computational methods developed within the 
theory of gases |5j| in a similar way to |^0| . 

In the instantaneous approximation we let particles to 
make a sudden jump from their in-coming trajectories 
into their out-going ones so that at time t + A t particles 
arrive at the correct coordinates, r' a and r' bl with the cor- 
rect momenta, k' and K' . Accordingly, in the asymptotic 
region, after t + A t , there is no distinction between the 
non-instantaneous and instantaneous pictures. 

The actual time instant at which the jump happens 
does not influence the asymptotic states, however, two 
particular choices of the time instant are important with 
respect to the implementation into simulations. The first 
one is the time at which the simulation code selects two 
particles as adepts for the collision usually at the point 
of closest approach. We will discuss this time later when 
we specify the simulation scheme. The second one is the 
central time, i = t + ^A t , for which the instantaneous 
approximation maintains the space and time symmetries 
of the non-instantaneous collision. As we will see, these 
symmetries allow one to derive the value of the sudden 
jump from experimental phase shifts. 

The condition that the sudden jump correctly mimics 
the non-instantaneous process is naturally met if one ex- 
trapolates the in-coming and out-going trajectories from 
known coordinates and momenta at t and t + At, re- 
spectively, to the central time t. Doing so one finds that 
extrapolated coordinates just before and after the sudden 
jump read 
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In extrapolations of momenta we neglect Coulomb forces 
so that we assume that protons and neutrons are driven 
by the same force from the Skyrme potential U, F = 
— ^f- In this case, the relative momentum remains un- 
changed and the force affects only the sum momentum, 



K = K + FA t , 
k = k = k — q, 
K' = K' - FA t = K 
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Finally, energies at the extrapolated phase-space points 
are 
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In rearrangement we have used that the increase of the 
kinetic energy, say £ - £ = ^ (fc + ±FA t ) 2 - £ = 
Y~ At, is compensated by the decrease of the poten- 



tial energy along the trajectory, U(f a ,t) — U(r a ,t) — 
|^2^At = "^^m^t- l n a stationary potential, the 
compensation of kinetic and potential energies reflects 
the energy conservation. The second order corrections in 
A's are neglected. 

Using the extrapolated quantities, we can define a new 
set of effective A's corresponding to the instantaneous 
picture, 
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The space displacements A bo , A a and A b can be ex- 
pressed in terms of the on-shell scattering phase shift 
defined as 
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From (H) and (||) one can directly check that 
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Effective displacements ( |ll[ ) can be evaluated from the 
experimentally observed scattering phase shifts. Before 
we turn to this pragmatic question, it is profitable to 
enlighten the conservation laws and symmetries of the 
collision processes. 



C. Conservation laws 

The extrapolated momentum and energy gains vanish, 
A K = K' - K = A K - 2FA t = 0, 
A E = ~e' a + e' b - 
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h = A E - 2— At = 0. (12) 



To show this we use that for a collision of isolated nucle- 
ons, the scattering phase shift depends only on the initial 
and final momenta, k and n' , while sum momentum K 
and the Skyrme potential only shifts the energy bottom, 

<t>(il, fc, p, q, r,t) = 4>(n- (fc +J )2 - 2U, k - p, q 

(13) 

The time derivative which results in the energy gain A e 
thus can be expressed via the energy derivative and from 
(|) one finds that A E = 2^A t , therefore A E = 0. The 
space dependency also enters (j) only via the energy argu- 
ment, therefore Ak = 2FA t or = 0. The cancella- 
tion of both shifts simplifies the energy and momentum 
conservation to its form commonly used for instantaneous 
collisions. Briefly, in the instant collision the Skyrme po- 
tential has no time to pass any energy and momentum 
to colliding nucleons. 

An additional simplification follows from the continu- 
ity of the centre of mass motion. This requires f' a + r' b = 
fa + fb or A a + Af, = 0. This relation is satisfies by dis- 
placemcntsfO) . Indeed, in the approximation of isolated 
collision (|l3|), the on-shell energy argument reduces, 
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The on-shell scattering phase shift ( |10| ) then does not de- 
pends on the sum momentum fc+p, therefore the deriva- 
tives with respect to fc and p are mutually connected, 
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According to (p. l|), the displacement of the particle b is 
opposite to the displacement of the particle a, 

A h = -A a . (16) 



D. Rotational symmetry 

The symmetries are best seen in the barycentric rep- 
resentation in which the phase shift is a function of the 
initial and final relative momenta, <p(k —p,q) = 4>(k, re'), 
where re — i(fc — p) and re' = ^(k — p — 2q). From sub- 
stitution into the barycentric framework one obtains 
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therefore the displacements in terms of relative momenta 
read 
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Apparently, the A a does not change under replacement 
of the initial and final momenta, n < — ► n' . As a comple- 
mentary displacement we have introduced A HS instead of 
A be which reverses its orientation under the replacement 
of initial and final momenta. 

For central forces, the scattering phase shift has to 
satisfy the rotational symmetry, therefore it depends only 
on the deflection angle 9 
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and amplitudes of initial and final relative momenta, 

4>(n,n') = 4>{cos 9, \ K \,\n'\). (20) 
The vector derivatives follow from 
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The derivative with respect to the final momentum n' is 



obtained from (|21l) via the interchange n 



E. Time-reversal symmetry 

As mentioned above, we want to express the effective 
displacements in terms of observable scattering phase 
shifts. This will make possible to circumvent the un- 
certainty about interaction potentials and to supply the 
simulation codes directly with experimental values. 

Although the phase shift </> does not depend on the 
energy O, it is still not experimentally known for general 
momentum q but only for momenta which satisfy the 
energy conservation, (k—p) 2 — (k—p—2q) 2 . Due to time- 
reversal symmetry of the collision process, it is possible to 
express the effective displacements in terms of derivatives 
along this momentum shell. 

The experimentally available values of the scattering 
phase shift are restricted to the shell \n'\ = \n\, 
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In (Eh) we need separate derivations with respect to \n\ 
and^K'|. Fortunately, due to the time and space rever- 
sal symmetries, the scattering phase shift is a symmetric 
function of |k| and \n'\, 
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Since we need to find the derivatives only for |re'| = |re| 
and symmetry (p3|) implies 
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we obtain a direct link to observable scattering phase 
shift 
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F. Displacements in simulations 

Now we are ready to evaluate the effective displace- 
ments from experimental scattering phase shifts. From 
(H), © and (H) we find 
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In agreement with their even/odd symmetry under in- 
terchange of initial and final momenta, the A Q is pro- 
portional to the sum re + re' while A HS to the difference 
re + re'. From the energy conservation, (re + re') (re — re') = 
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k 2 — k' 2 = 0, follows that these two vectors are orthogo- 
nal. Before we present actual values of displacements, it 
is profitable to compare their form with the model of hard 
spheres and the approximation by the collision delay. 

For classical hard spheres of radius R, the scattering 
phase shift depends only on the transferred momentum 
q = k — k' as, </> HS = 7r — 2|g|i?. From ( prf ) than follows 
that A a = and A HS = ^2R. The hard-sphere approx- 
imation used in model studies Q [3(| thus neglects A a 
and uses a constant approximation of the amplitude of 



the other displacement, 



= 2R. This amplitude is 



conveniently evaluated from the cross section, a = irR 2 
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FIG. 2. The effective displacement as a function of the de- 
flection angle and the kinetic energy, K 2 /m, in the barycentric 
coordinate system. The columns of dots show the spread of 
components with deflection angle. The lines show the an- 
gle-averaged values. The amplitude of the T-matrix is pre- 
sented in the top section to indicate the weight of individual 
processes. The orthogonal component, Ax/2 = d(f>/dcos8/2K 
shown in the middle section, has appreciably smaller values 
than the parallel component, Ay /2m = d<f>/dn/2 shown in 
the bottom section. 

The approximation discussed by Danielewicz and Pratt 
f33f deals only with the collision delay defined according 
to Wigner as the energy derivative of the scattering phase 
shift. This approximation is obtained from ( p6| ) if one 
neglects the derivative with respect to the deflection an- 
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where E = \k\ 2 /m is energy in the barycentric system. 
Numerical values of these two contributions are com- 



pared in Fig. The dots in the vertical line show 
a spread of values due to the angular dependence, the 
curves show values averaged over deflection angles with 
the weight given by the differential cross section displayed 
in the top section. The parallel component, shown in 
the bottom section, has a typical value of 0.5 fm. The 
negative large values below 3 MeV can be ignored since 
corresponding processes have very small rates due to the 
Pauli blocking. The perpendicular component, shown in 
the middle section, has about three-times smaller values, 
moreover it tends to average out. For energies above 
lOMeV, the displacements can be well approximated by 
a constant value, as it is case of the hard-sphere model. 
Moreover, the amplitude of the displacement is close to 
the estimate based on the differential cross section, in 
spate conceptual difference between both concepts. Our 
results thus confirm that estimates used in [|^|- p0|j5lj| are 
quite reasonable. 

In f?8[ we have used as the time of instant jump the 
time of closest approach. This distance is different from 
the distance A bc required from the equivalent scattering 
scenario presented in figure |^ as solid line. We consider 
now the time required to travel from Ab e to the dis- 



tance of closest approach A t 
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|36[ . Within this scenario we are allowed to jump at the 
point of closest approach to the final asymptotics ^ and 
(j?])with the additional distance the particle travel during 
At. The effective final jump in the barycentric frame is 
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where the same compensation of off-shell derivatives by 
At occurs as described before when jumping at the centre 
time t. Please remark that (ETT) agrees with (§). 
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FIG. 3. A nonlocal binary collision (solid line) together 
with the scenario of sudden jump at the closest approach. 



G 



G. Renormalisation of quasiparticle energies 

So far we have discussed the nonlocal shifts as if there 
were free classical particles. The interaction affects, how- 
ever, the free motion of particles between individual col- 
lisions. The dominant effect is due to mean-field forces 
which bind the nucleus together, accelerating particles 
close to the surface towards the centre. These forces 
are conveniently included via potentials of Skyrme and 
Hartree type. 

Beside forces, the interaction also modifies the velocity 
with which a particle of a given momentum propagates in 
the system. This effect is known as the mass renormali- 
sation. A numerical implementation of the renormalised 
mass is rather involved since a plain use of the renor- 
malised mass instead of the free one leads to incorrect 
currents. Within the Landau concept of quasiparticlcs, 
this problem is cured by the back flow, but it is not ob- 
vious how to implement the back flow within the BUU 
simulation scheme. In our studies, we circumvent the 
problem of back flows using explicit zero-angle collisions 
to which we add a non-local correction. 



1. One- dimensional system of fixed scatterers 

A link between the renormalisation of the mass and 
zero-angle collision has been already pointed out by Lan- 
dau. We find it instructive to describe this mechanism 
first for a simple one-dimensional system of randomly 
distributed barriers. Tunnelling then corresponds to the 
zero-angle scattering and reflection to a dissipative colli- 
sion. We focus on the piece of quasi- free trajectory, i.e., 
the trajectory between two successive dissipative colli- 
sions. 

A tunnelling through barriers speeds up or slow down 
the mean velocity of particles. To simulate this effect 
on the motion, at each tunnelling we shift a particle by 
a displacement D in the direction of its motion. For 
simplicity we take the amplitude of D as constant. After 
time t the particle moves over a distance 



k 

x=t \-ND, 

m 



(28) 



where k/m is its velocity between tunncllings and N is 
the number of barriers on the trajectory of length x. In 
average N = cx, where c is concentration of barriers. 
The mean velocity of a particle then reads 



x k „ x k 

v = - = h cD- = — . - 

t m t m(l — cD) m 



— {1 + cD). (29) 



For the three-dimensional system with particle-particle 
interactions, the corrections to the mean velocity will be 
limited to the linear approximation. 

Relation (|2^) can be compared with the velocity of a 
particle evaluated within the renormalised mass m*, 



in" 



from which follows 



1 + cD. 



(30) 



(31) 



We can use (^lj) to fit D so that the known renormalised 
mass m* is reproduced. 



2. Three-dimensional Fermi 



In the classical three-dimensional system all collisions 
have a finite deflection angle. In the quantum system, 
however, there are zero-angle collisions which represent 
an interference between scattering states and the incom- 
ing state of the interaction. In the dialect of the pertur- 
bative expansion one can say that the particle makes a 
detour from its trajectory in the phase space but nowhere 
on the detour it reaches the energy shell and thus it has 
to return back. The detour causes a delay expressed by 
the shift D. In the Fermi liquid, the Pauli exclusion prin- 
ciple blocks a majority of phase space cell on the energy 
shell so that the zero-angle collisions dominate over the 
dissipative events. We will use these blocked events to 
simulate the renormalisation of the mass. 

Unlike in the simple one-dimensional scattering on 
fixed defects, the displacement D is a vector oriented 
along the difference momentum, 



D = \D\ 



P 



\k 



■P\ 



(32) 



where k is a momentum of the assumed particle while 
p belongs to its partner in the prohibited collision. In 
general, the displacement \D\ is a function of k and p. 
For simplicity we assume this function as a constant and 
fit its value to the mass renormalisation at the Fermi 
surface. 

For the fitting of the displacement D we assume the 
zero temperature at which all real collisions are blocked 
by the Pauli exclusion principle so that all binary en- 
counters contribute to the renormalisation. The mean 
velocity of the particle is then given by its free motion 
and the mean value of the displacements per the time 
unit, 



k 

v = — 
m 



E 



dp \k-p\ k 

a fp\ D \Tr 



p 



in 



■P 



(33) 



The mean value of displacements is proportional to the 
frequency of binary entertainments, i.e., it is the sum 
of integrals over distributions of protons and neutrons 
weighted with the scattering cross section a and their 
relative velocity to the observed particle. 

With a good approximation the cross section a is in- 
dependent of energy so that one can easily evaluate the 
integral in (p3|), 
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k k 

v = h na\D\ na\D\(v). 

m ui 



(34) 



The renormalised velocity thus depends on the density, 

(35) 



3w 2 h 3 ' 

and the mean velocity of the nuclear matter 



n J (^Trn) 6 m 



P 



(36) 



For the system in rest, (v) = 0, we find the mass renor- 
malisation, 



= l + na\D\ 



(37) 



One can see that this formula has the same interpreta- 
tion as the one-dimensional case (^lj), because na is the 
average number of scatterers on the trajectory of unitary 
length. 

Formula ( |37| ) allows us to fit |D| from known value of 
the effective mass. For m* : m — 3 : 4, a — 40 mb 
and n = 0.16 fm -3 one finds value \D\ = 0.5 fm. This 
value is very close to the nonlocal correction in dissipative 
collisions, see Fig. ||. 

The quasiparticle velocity (|34|) relates to the quasipar- 
ticle energy by v = dek/dk. For the moving nuclear 
matter (v) ^ one finds 



k 2 (k — m(v)) 2 

e k = -z \-na\D\ . 

2m 2m 



(38) 



An approximation of this structure is commonly used in 
simple applications of the Landau concept of quasiparti- 
cles. 



not least, the simulation does not require to introduce 
new time-demanding procedures, one can simply use the 
scattering events which are merely rejected in standard 
simulation codes by Pauli-blocking. 



III. NUMERICAL RESULTS 

Let us discuss the proposed correction to the local 
and ideal (no quasiparticle renormalisation) Boltzmann 
(BUU) simulation. First we introduce the pure nonlocal 
corrections and then we discuss the quasiparticle renor- 
malisation. 

The evolution of the density can be seen in the corre- 
sponding left pictures of figure | for the BUU (left panel) 
and nonlocal scenario (middle panel) as well as the ad- 
ditional quasiparticle renormalisation (right panel). We 
see that the nonlocal scenario leads to a longer and more 
pronounced neck formation between 200 — 240fm/c while 
the BUU breaks apart already at 200fm/c. 

The question arises whether this pronounced neck for- 
mation is simply by more collisions and corresponding 
correlations. This would lead us to the assumption that 
a simply increase in the cross section as sometimes called 
in-medium effect would lead to the mid-rapidity matter 
enhancement. This is however only the case for smaller 
impact parameters pi| . To understand the qualitative 
difference between the nonlocal scenario and an increase 
of cross section by in-medium effects wc perform a simu- 
lation where in the local BUU scenario the cross section 
has been doubled. We see in in the next figure || the 
number of collisions per time for the different scenarios. 



H. Summary an simulation schema 

The derived effective nonlocal collision procedure is 
easily incorporated in the usual collision simulation by 
an additional advection step. The quasiparticle renor- 
malisation and effective mass is found to be possible to 
incorporate by the same advection step but performed 
for the events rejected normally by Pauli-blocking. 

Finally, we would like to comment on properties of 
the proposed simulation scheme. The renormalisation 
depends on the distribution of particles in surrounding 
medium. It has four nice properties: (i) the renormali- 
sation vanishes as the local density goes to zero, (ii) the 
renormalisation vanishes when a high temperature closes 
the Luttinger gap because all collisions will be at finite 
angles, (iii) the anisotropy of the quasiparticle velocity 
in a presence of a non-zero current in medium is auto- 
matically covered, and (iv) the backflows connected to 
the mass renormalisation are covered because both par- 
ticles jump keeping the centre of mass fixed. Last but 
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FIG. 4. The time evolution of \f-Ta +\% 7 Au collisions at Ei ab /A = 33 MeV and 8fm impact parameter in the BUU (left), 
nonlocal kinetic model (middle) as well as the nonlocal model with quasiparticle renormalisations (right). Plots in the first 
column show the (x — z) - density cut where Ta as projectile comes from below. The mass momenta are shown by arrows. 
The corresponding second column gives the charge density distribution versus relative velocity in cm/ns where the target like 
distribution of Au is on the left and the projectile like distributions of Ta on the right. 



10 





40 




30 




20 


fm] 


1 n 

I u 


^5 





o 




E 


30 






o 


20 


Q. 




CO 
C 


10 


o 




« 





o 


30 


o 




20 




10 








A b=0 fm 






b=6 fm 






b=8 fm 


b=12fm 







u(r, i) 



1 



100 



200 



100 



200 



300 





t [fm/c] 

FIG. 5. The time evolution of the number of nucleon col- 
lisions for Ta + Au at Ei a b/A — 33 MeV and different im- 
pact parameter in the BUU (thick black line), nonlocal kinetic 
model (broken line) and for the case of 8fm impact parameter 
the local BUU with a cross section twice as large (thin dark 
line) . 

For 8 fm impact parameter we compare the local BUU 
(thick black line) with the nonlocal (broken line) and the 
local BUU multiplying the cross section with two (thin 
line). We see that the number of collisions are visibly en- 
hanced by doubling the cross section while for the nonlo- 
cal scenario we get only a slight enhancement at the be- 
ginning and later even lower values with respect to local 
BUU. The latter fact comes from the earlier decompo- 
sition of matter in the nonlocal scenario. Consequently 
from the number of collisions we would conclude that the 
increase of cross section leads to more correlations than 
the nonlocal scenario. 

However when we look at the corresponding transverse 
and kinetic energies in figure || (8fm impact parameter) 
we see that the transverse and longitudinal energy is al- 
most not changed compared with local BUU. Oppositely 
the nonlocal scenario leads to an increase of transverse 
energy of about 2MeV and about lMeV in longitudinal 
energy. We conclude that the increase of cross section 
leads to a higher number of collisions but not to more 
dissipated energy while the nonlocal scenario does not 
change the number of collisions much but the energy dis- 
sipated during the collisions. Roughly speaking we can 
say that the quality of collision is changed. 

Returning to the discussion of pronounced neck forma- 
tion in figure |i| above we see now that the quality rather 
than the quantity of collisions is what produces the neck. 
The simple increase of the number of collision does not 
change much. 

Now we can proceed and discuss the charge matter 
distribution with respect to the velocity. We evaluate 
the mean density and velocity, 



dp 

mn(r, t) J (27r) ; 



pf(p, r, t), 



(39) 



from which we define the distribution of hydrodynamical 
velocities 



F(v,t) — / drn(r,t) 5(v - vti ss (r,t)) 



(40) 



where va ss (r, t) is the projection of v{r,t) onto the fis- 
sion line. This distribution we identify with the so called 
charge density distribution. 

The definition of mean mass (current) velocity does 
not include the Fermi energy which is integrated out. 
In the case that we do have a different repartitioning of 
Fermi energy during the collision than described in our 
kinetic equation we will have here an ambiguity. Since 
the dynamical cluster formation is not described in our 
approach we might have here a smaller effect of Fermi 
energy on the mass velocity. This will lead us indeed 
to the observation that BUU or nonlocal kinetic equa- 
tions have too much stopping compared to the experi- 
ment when more central collisions are considered. For 
peripheral collisions we believe that this kinetic descrip- 
tion is sufficient which we will prove by proper association 
of experimental events to the maximum in the velocity 
distribution. 
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FIG. 6. The time evolution of the longitudinal (thin lines) 
and transverse energy (thick lines) including Fermi motion 
of nucleon collisions for Ta + Au at Ei a b/A = 33 MeV and 
different impact parameter in the BUU (black line), nonlocal 
kinetic equation (dashed line) and for the case of 8fm impact 
parameter the local BUU but with twice the cross section 
(dotted line). 

We plot in figure [| also the normalised charge dis- 
tribution versus velocity and see that after 160fm/c we 
have an appreciable higher mid-rapidity distribution for 
the nonlocal scenario (mid panel) than the BUU (left 
panel). Together with the observation that for nonlocal 
scenario we have a pronounced neck formation we see in- 
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deed that the neck formation is accompanied with high 
mid-velocity distribution of matter. 

A. Quasiparticle renormalisation 



This softening of equation of state is already slightly re- 
markable when the nonlocal scenario is compared with 
BUU. With additional quasiparticle renormalisation we 
see that this is much pronounced. 



Now we use the quasiparticle ren orma lisation schema 
which has been outlined in chaptei HG . We see in fig- 
ure ^ (right panel) that the mid-rapidity distribution of 
matter is once more enhanced in comparison to nonlo- 
cal scenario without quasiparticle renormalisation. The 
seemingly shorter lifetime of the neck is artificial due to 
the chosen density contours which means we have lower 
densities and faster matter disintegration so that in fact 
the neck is much more pronounced than in simple non- 
local scenario and of course more pronounced than in 
BUU. The detailed comparison of the time evolutions of 
the transverse energy for 8fm impact parameter can be 
seen in figure 0. 
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FIG. 7. The time evolution of the transverse energy includ- 
ing Fermi motion for Ta + Au at Ei a h/A = 33 MeV and 8fm 
impact parameter in the BUU (black line), nonlocal kinetic 
equation (dashed line), the local BUU but twice cross section 
(dashed dotted line) and the nonlocal scenario with quasipar- 
ticle renormalisation (long dashed line). 

We recognise that the transverse energies including 
quasiparticle renormalisation are similar to the nonlocal 
scenario and higher than the BUU or BUU with twice 
the cross section. However please remark that the period 
of oscillation in the transverse energy which corresponds 
to a giant resonance becomes larger for the case with 
quasiparticle renormalisation. Since therefore the en- 
ergy of this resonance decreases we can conclude that the 
compressibility has been decreased by the quasiparticle 
renormalisation. Sometimes this quasiparticle renormal- 
isation have been introduced by momentum dependent 
mean-fields. The effect is known to soften the equation 
of state. We see here that we get a dynamical quasiparti- 
cle renormalisation and a softening of equation of state. 



B. Comparison with experiments 

The BUU simulations will now be compared to one ex- 
periment performed with INDRA at GANIL, the Ta+Au 
collision at Eiab/A — 33 MeV. The first question when 
comparing with experiments concerns the proper selec- 
tion of events such that one can compare with specific 
impact parameter of the simulation. We choose here the 
point of view that the maximum in the charge distri- 
bution with respect to velocity which is a measure for 
stopping gives a good correlation with impact parame- 
ter. Indeed if we compare the corresponding correlation 
between impact parameter and this maximum velocity 
we obtain indeed an almost linear correlation as in figure 
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FIG. 8. The maximum velocity, impact parameter and ra- 
tio of longitudinal to total kinetic energy of Ta + Au collisions 
at Eiab/A = 33 MeV in the BUU (solid line) and the nonlocal 
model with quasiparticle renormalisation (dotted line). 

The matter distribution is shown for different approxi- 
mations in figure g. One recognises clearly the successive 
enhancement of mid-rapidity matter around 6 — 8fm if 
one uses nonlocal kinetic theory and quasiparticle renor- 
malisation correspondingly. It is interesting to remark 
that the dynamical quasiparticle renormalisation which 
leads to a softening of the equation of state as discussed 
in figure enhances the mid rapidity distribution. In con- 
trast a mere soft static paramctrisation of the mean-field 
does not change the mid-rapidity emission appreciably 

For the identification with experimental selection we 
use the selection of events in the following way. First we 
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select events which show a clear one fragment structure. 
This correspond to events where we have clear target and 
projectile like residues. Since the used kinetic theory is 
not capable to describe dynamical fragment formation we 
believe that these events are the one which are at least 
describable within our frame. Next we use impact param- 
eter cuts with respect to the transverse energy since this 
shows in all simulation a fairly good correlation to the im- 
pact parameter. In our numerical results we see almost 
linear correlations between impact parameter, maximal 
velocity and the convenient ratio between transverse and 
total kinetic energy as seen in figure |^. 

For each selected experimental transverse energy bin 
we can plot now the maximum velocity versus the ra- 
tio of the transverse to kinetic energy. We see in figure 
[To| that the numerical velocity damping agrees with the 
experimental selection only for very peripheral collisions. 
For such events we plot in the figure [l^ the charge density 
distribution and compare the experiment with the sim- 
ulation. These charge density distributions have been 
obtained using the procedure described in reference J45| . 
The Data are represented by light grey points, the stan- 
dard BUU calculation by the thin line and the non-local 
BUU with quasiparticle renormalisation calculation by 
the thick line. A reasonable agreement is found for the 
nonlocal scenario including quasiparticle renormalisation 
while simple BUU fails to reproduce mid-rapidity mat- 
ter. 
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line) as well as the nonlocal model with quasiparticle renor- 
malisation (dotted line). 



compared to recent experimental data on mid rapidity 
charge distributions. It is found that both the nonlocal 
shifts as well as the quasiparticle renormalisation must 
be included in order to get the observed mid-rapidity 
matter enhancement. 

The inclusion of quasiparticle renormalisation has been 
performed by using the normally excluded events by 
Pauli blocking. Since the quasiparticle renormalisation 
and corresponding effective mass features can be consid- 
ered as zero angle collisions they can be realized by non- 
local shifts for the scattering events which are normally 
rejected. This means that one has to perform the advec- 
tion step for the cases of Pauli blocked collisions without 
colliding the particles. Besides giving a better descrip- 
tion of experiments, this has the effect of a dynamically 
softening of equation of state seen in longer oscillations 
of giant compressional resonance. 

In this way we present a combined picture including 
nonlocal off-sets representing the nonlocal character of 
scattering, which leads to virial correlations with the 
quasiparticle renormalisation, and as a result to mean 
field fluctuations. We propose that no additional stochas- 
ticity need to be assumed in order to get realistic fluctu- 
ations. 
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IV. SUMMARY AND CONCLUSION 



The extension of BUU simulations by nonlocal shifts 
and quasiparticle renormalisation has been presented and 
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FIG. 10. The experimental charge distribution of matter (dotted line) versus velocity in comparison with in the BUU (thin 
solid line) and the nonlocal model with quasiparticle renormalisation (thick line). The maximum velocity versus ratio of 
longitudinal to total kinetic energy of Ta + Au collisions at Ei a b/A = 33 MeV is given below. The selected experimental cuts 
are given by dots. 
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